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TECHNICAL MEMORANDUM 


MINIMUM WEIGHT DESIGN OF A LEAF SPRING TAPERED IN THICKNESS AND 
WIDTH FOR THE HUBBLE SPACE TELESCOPE-SPACE SUPPORT EQUIPMENT 

INTRODUCTION 


During the life of the Hubble Space Telescope (HST). on-board optical guidance systems 
and scientific instruments will experience degradation. Maintenance or replacement of these systems 
will be necessary in order to maintain a fully operational observatory. Due to the cost and risk of 
retrieving the HST for ground refurbishment and consequent space redeployment, a series of main- 
tenance missions have been identified in order to carry fine guidance sensors (FGS's) and scientific 
instruments (Si's) aboard the space shuttle for on-orbit replacement of degraded units. The weight 
of these units ranges from approximately 500 to 1.000 lb. 

During the initial launch of the HST. these instruments form part of a 25 .000-1 b space 
observatory. This large mass provides the Si s and FGS's with a safe environment from the fre- 
quency spectrum of the space shuttle cargo bay. When launched separately as part of a mainte- 
nance mission, the protection from the dynamic environment must come from a suspension system 
that will preclude damage to these delicate optical and scientific instruments. 

As part of the design of the suspension system, leaf springs (similar to those found in 
automobiles) have been designed to provide the necessary flexibility to alleviate potentially damag- 
ing dynamic loading. This report describes the design of a concept of a variable width and depth 
cantilever spring for the HST maintenance and refurbishment mission. 


LINEAR ELASTIC SOLUTION 


The basic idea behind the suspension system is to provide a certain stiffness (or flexibility) 
so that there is no danger of resonance between the natural frequency of the system (including the 
payload) and certain mechanical and acoustical frequencies encountered during the ascent or descent 
phase of the mission. There are also certain points in the ascent/descent frequency spectrum that 
could induce high transient loads into the hardware and. thus, the natural frequency of the suspen- 
sion system should be different than these. 

The first step in the design of the spring is to determine its stiffness or spring rate. Because 
the design of the suspension system will limit the movement of the spring/mass system to one 
direction, its natural frequency can be obtained from the equation of motion of a single degree-of- 
freedom system 


ni\ + k.x = 0 


( 1 ) 


Solution of equation (1) leads to definition of the natural frequency of the system [1], 


fn = 


2it 



( 2 ) 


and solving for the stiffness k one obtains 


k = 4m(/ir) 2 . (3) 

With the stiffness k and the mass m to be isolated, one has the information necessary to begin the 
design of the leaf spring. In order to produce an efficient spring that is compatible with the avail- 
able installation space, a cantilever with variable width and depth has been selected. This is a 
hybrid between type F-3 (triangular cantilever) and type T-l (tapered cantilever) which can be 
found in chapter 10 of reference 2. Figures 1 and 2 show the suspension system assembly and the 
geometry of the spring, respectively. At this point, it should be explained that the beam in figure 2 
has the load applied at the two tapered ends through pins. The center (constant cross section) of 
the beam is clamped and bolted to essentially provide two cantilever springs instead of a longer 
simply supported spring. 

By defining the widths and thicknesses at the fixed and free ends of the spring as w„, \v e 
and t (> , t ( „ respectively, one can express a linear taper for both the width and depth as 


(f) +, '"( 1- f) • 

(4) 

(t) + '"( ,_ r) • 

(5) 


where tv(jc) and t(x) are the width and thickness at any position .v along the length L of each canti- 
lever (the origin x„ is located at the fixed (clamped) end). The moment of inertia at any point 
along the length is 


/(.v) = 


ntvm.r) 3 

12 


(6) 


Substituting equations (4) and (5) into equation (6) one obtains 


/(A) = 


(U’|.V + H’tH/ |.Y + f t )' 


12 


(7) 


where 


“’I 


u\. - u\ 


L 


( 8 ) 


2 




3 


( 9 ) 


VV’2 = U', 


'l 



h = t, 


( 10 ) 

(ID 


The relationship between the elastic axis of the spring and the bending moment for elastic 
deformations is given by the Euler-Bernoulli equation. 


J_ _ Mix) 
P EKx) 


(12) 


Equation (12), although derived for prismatic bars, can be utilized in the analysis of tapered beams 
with sufficient accuracy as long as the variation of the taper is not extreme f3]. The equation for 
the curvature of the elastic beam is 


j_ = (l 2 y/ d\ 2 

0 1 1 + Uty/dx) 2 } 212 


(13) 


When dealing with small deflections, which correspond to the linear range of equation (13), the ef- 
fect of the dy/d.x term becomes negligible and one can write 

(1 + (dy/d.x) 2 ] 212 « 1 . (14) 

Combining equations (12), (13), and (14), one obtains the linear elastic range of the Euler- 
Bernoulli equation. 


d 2 y _ M(-v) 
dx 2 El(.x) 


(15) 


Referring to figure 4 and recalling that the moment of inertia is a function of .v. one can 
express equation (15) as 


n d 2 x _ 12 P(L-.x) 

b — r 

dx~ (u’|.v + u’ 2 )(/|.v + /t)' 


(16) 


where 


M(.x) = P(L-.x) 


(17) 
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and /(.v) is defined as in equation (7). The right hand side of equation (16) can be expressed as a 
sum of partial fractions as follows. 


\2P(L-.x) _ A , B , C f D ( i8) 

(u |.v -Mr : )l/|.V + M' »V|A'+tV 2 O-'+M U\X + t 2 ) 2 (/|.V + / 2 ) 3 

The following are defined. 


a = 

(19a) 

b = 3/| 2 / 2 

(19b) 

v = .W/, 

(19c) 

(1 = /y 

( I9d) 

c = r t 2 ir. 

( I9e) 

ir 2 /| 2 + 2 /i/ 2 u - | 

(190 

Mi’i + 2 1 1 r 2 ir 2 

(I9g) 

/; = / 2 2 us 

( 19h) 

j = 0"'i 

(I9i) 

= 1 1 us + r 2 u'| 

( 1 9 j ) 

= Mis 

(19k) 


By solving the partial fractions problem of equation (18). and by using equations (19). one can 
obtain the expressions for the constants A. B. C, and D. In matrix form, these equations are. 



( 20 ) 
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Substituting the values of A, B, C, and D obtained in equations (20) into equation (18), one can 
now express equation (16) as follows, 


E^ = + - + — - + - . (21) 

dx~ h’|.v + vi '2 r ,.y + 1 2 V |.v + / 2 ) - (ri.v + fi) 3 

Performing consecutive integrations on equation (21) will lead to the equation for the deflection of 
the beam. The first integration leads to the equation of the slope of the beam. 

E f ( c -^r\ dx = E c p- = — ln(iV|.v + «’■») + — ln(/|.v + f->) 

L \ d\~ ) dx h'i - 1 1 

( 22 ) 

- — - — - + d> . 

/i(r|.v + r 2 ) 2ri(/|.v+/2)" 

<f> is a constant of integration, which evaluated at .v = 0 with dy/dx = 0 yields 


<t> = -£-+■ Q. -A \ nW y—lL in/ 2 . (23) 

hh 2 r l t2~ u ’i 

Performing a second integration on equation (22) leads to the expression for the deflection at any 
point on the beam. 


E {(£) dr = Ey = >^[''' ! ^ {ln<H ' | - t+u - )>_ ' v ] + f‘ +'=»-*] 

.£[i w , It „)ll + ff^ — -1 
h L h J 2t\ / 1 (/ |.v + / 2 ) J 


+ <J>.v + T , 


(24) 


w 


here V P is a second constant of integration which can be evaluated at v = 0 where y — 0. 


— — — — (itslnns) — — -(r 2 lnt 2 ) T ~ ~( ln/ 2 ) — — — (25) 

n-|- 1 1" t\~ 2f|-r 2 

Equation (24) is the linear elastic solution for the deflection of a variable thickness and width 
cantilever beam. 

The maximum bending stress for the beam of constant cross section will occur at the loca- 
tion where the bending moment is maximum. This is not the case for a beam having a variable 
thickness. Since the moment of inertia of the beam is a function of the cube of the thickness, the 
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section modulus changes much faster than the bending moment as a function of beam length. This 
causes the maximum bending stress to occur at a location other than the point of maximum bend- 
ing moment. A normalized plot comparing the bending moment, bending stress, and moment of 
inertia for a typical cantilevered beam of tapered width and thickness can be found in figure 3. The 
expression for bending stress for the beam under study is 




6 P(L-.v) 

(u’|.v + \is)(r |.r + /t)~ 


(26) 


Since the maximum bending stress is now a function of the varying cross section as well as 
the location along the length of the beam, one can find the maximum by setting the derivative of 
equation (26) equal to zero. 


= 0 . (27) 

ax 


The resulting expression is a cubic equation in x which can be solved to obtain the location ol 
maximum stress 


where 


.v 3 + a ,.v 2 + an.v + a-, = 0 , 

(28) 


(29a) 

“ 3 = - i [ 2 0r) + *d • 

(29b) 

T ns t->L . t-rL . \\’ 2 h 2 

|_ w’i'i 2t, 2 2 m -,/, 2 J 

(29c) 


All three roots of equation (28) will be real, with only one root being physically meaningful. 
Substitution of the location of maximum stress ,v max (feasible root of equation (28)) into equation 
(26) yields the maximum stress in the beam. At this point, one has the necessary information to 
design a beam with a desired stiffness and subject to a maximum allowable bending stress. 

WEIGHT OPTIMIZATION 


In the aerospace industry, minimum weight of flight structures is of primary importance to 
the structural design engineer. The use of the space shuttle as a space transportation vehicle results 
in an approximate cost of $1,100.00 per pound to deliver a payload to Earth orbit. It is. therefore, 
obvious that lighter payloads result in lower costs and are in the best interest of the government. 
The challenge is to be able to minimize the weight of a design by varying a select group of 
parameters and not violate constraints that are essential to the structural integrity of the hardware. 
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Normalized Variables 


00 



Figure 3. Normalized bending stress variation lor a typical double tapered cantilever beam. 


During recent years, advancements in the field of mathematical programming coupled with 
the everchanging state-of-the-art in personal computer hardware have provided a great opportunity 
for improvement in optimization software and code availability. Commercial and academic optimi- 
zation packages are currently available that will work on almost every personal computer hardware 
platform. Examples of these are I-DESIGN (University of Iowa. Dr. Jasbir Arora) and DOT 
( Vanderplaats. Miura & Associates). Other packages that have been generated under government 
contracts are ADS. CONMIN. ACCESS, and NEWSUMT. Additional information on optimization 
software availability can be found in references 4. 5. and 6. 

The general problem statement for the minimization of a function of several variables 
subject to conditions of constraint is. 


Minimize: 

F(X) 

(30) 

Subject to: 

g/(A) 0 

j = 1 - nt 

(31) 


/t*( X) = 0 

k=\.l . 

(32) 


where X is the vector containing the design variables. F is the objective function (function to be 
minimized), g, are the inequality constraints, and /; k are the equality constraints. In order to limit 
the region of search for the optimum, side constraints are imposed on the problem. This is accom- 
plished by simply imposing upper and lower bounds on the search values of the design variables. 

X/^Xi^Xi" /'= 1./? . (33) 

For the problem of weight minimization the objective function is the volume of the beam. This 
volume can be obtained by integrating a differential element of area over the entire length of the 
beam. 


Vol = f u i.v)t(.x)dx . (34) 

0 

where 


U’(.v) = U’j.V+U ’2 

(35) 

t(.\) = /|.V /2 

(36) 


and u - |. us. 1 1 . and i 2 are defined in equations (8) through (II). After performing the necessary 
integration, the objective function can be expressed as 

W(X) = A |.v( I ).v(2) + B | |.v( 1 ).v(4) + .v(2).v(3)] + C ,.v(2).v(4) . (37) 
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where 


= pL-V.3 

(38) 

= pL : /2 

(39) 

= p L . 

(40) 


and p is the density of the material. In order to express the design variables in a logical and con- 
Mstent manner for computer implementation, they have been identified as follows. 


.v( 1 ) = ll| 

(41a) 

v(2) - us 

(41b) 

.v( 3 ) = /, 

(41c) 

v(4) = t 2 

(4 1 d ) 


For the purpose of this report, it is desired to design a beam of minimum weight that does 
not exceed the allowable y ield stress of the material. This means that the working stress, equation 
1 20). ma\ not exceed the yield stress of the material. One can identify this restriction as an 
inequality constraint and. in normalized fashion, it can be expressed as follows. 


<T> 

o'. 


=£ 0 


(42) 


Since one is looking for a specific stiffness of the beam, the maximum deflection must be set 
equal to a prescribed value. This value determines the desired natural frequency of the beam for a 
prescribed load. Using equation (24). one can identify this restriction as an equality constraint and. 
normalized, it will be expressed as follows. 


- 1=0. (43) 

ball 

In equation (24). y indicates the deflection of the cantilevered beam at any location .v along 
its length. In equation (43). 8 is the maximum deflection of the beam which occurs at the tip of 
the cantilever. 8 a n is the deflection associated with the desired natural frequency of the system to 
be dynamically isolated. 

Minimization of equation (37). subject to constraint equations (42) and (43). identifies the 
optimization problem. This problem can be stated as follows: "Find the minimum weight W of a 
variable cross section beam under a specific loading condition that will not exceed the allowable 
material yield stress oy and will have a maximum deflection of 8 ;lM ." 
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The optimization software used in the solution of this problem is Design Optimization fools 
(DOT's). Version 2.00 of this commercially available software allows solution of the problem 
using two known methods, “modified method of feasible directions" (MMFD) or “sequential linear 
programming" (SLP). The MMFD is a modification of the method of feasible directions (MFD) in 
which equality constraints can be handled by including them as part of a pseudo-objective function. 
The MFD algorithm is not capable of effectively dealing with equality constraints. Description of 
the MFD and SLP algorithms can be found in references 5 and 7. Description of the MMFD can 
be found in reference 7. 


COMPUTER IMPLEMENTATION 


The cantilever beam of variable width and depth with deflection and stress constraints can 
nave many feasible solutions. Physically, this means that many combinations of »•„. and I, 

will lead to improved designs. Numerically, this means that the initial design variables must be 
carefully selected. The fact that the problem has many relative minima indicates that small vari- 
ations in the initial choices of design variables can lead to significant improvement in the design or 
nonconvergence. 

In an effort to provide reasonable first choices of initial variables, a computer program has 
been generated that will provide feasible design solutions. These solutions, although not necessarily 
the least weight designs, provide design variables that will attempt to meet the stress and deflection 
constraints. The program SPTRIAL thus provides initial solutions to the optimization problem. 

Once the initial values of the design variables have been chosen, they are input into the 
computer program SPOPT (SPring OPTimization) to obtain a design of minimum weight. SPOPT 
is a calling program that accesses the DOT optimizing software. Tables I and 2 are listings of the 
Fortran programs SPOPT and SPTRIAL. respectively. 

NUMERICAL EXAMPLES 


The initial design of the springs for the HST/space support equipment (HST/SSE) was 
performed without the benefit of optimization software. This means that numerous hand calcu- 
lations were performed and small, tailored computer programs were developed to aid in the many 
iterative calculations involved. The goal here was to obtain a design that would meet the stress and 
deflection constraints imposed. Weight minimization, although a big driver, was aimed at changing 
material/spring configuration combinations and not at refining the final geometry. Vast 
improvements in weight were accomplished by changing from coiled to multileaf to single leaf 
springs during the preliminary design phase. Once a configuration was selected, the refinement was 
limited to adjusting tolerances to meet the desired deflection while maintaining the stresses under 
the allowables. 

In this section, the author will start from the final spring geometry that resulted from the 
preliminary design phase and attempt to optimize the weight by using the SPTRIAL/SPOPT 
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Table I . Fortran program SPOPT. 


mien mn pricisior in mi 
di/hnsion mum, 5),xi wjuw, gw, aim 
amm mi m. iKiw.mm. mmi. smso/ 
mu mu. mi = X: i immooi dad 
opin(unii=b.fiii= 1: \mm\m bad 
m mn mi-- X: mmwmiMD 
c 

mm 
mmm 
DOWN. 10 
mu D=oo 

10 mm 
imu 
irihoch 

HUM 

NC0N=S 

C 

I/R 1 H 16 , •) xm m Blurt min' 
riadis, mo 
xm 

ipiii(6. •) mm wmi guiss for m mrns.m ■ 

Rim. •nil 

mu i6, •) mm Mim. guiss for no mimio.oi • 

Rim, •mo 

mn is, •/ mm mum guiss for n imimi m ■ 
Rim, •in 

mu (6, •) mm mum guiss for ;o mm= i. 50i • 
Rim, •no 

C ••••• MUM mils OF IHl 0FS1GR WIABLIS ••••• 

XU MU 

X/SMO 

X/JJ-IF 

X14U0 

C ••••• l (MR BOUNDS ON IHl 0FS16N XARIABKS ••••• 

XI It 1=999 
XI 121- 1.00 

XI (5 I =.2499 
XI (4) =15 

C ••••• IPPIfi BOUNDS ON 1HI DISIGN VAR1MHIS 
Ml 1=9. 999 

XU12/--I0.0 
XU(S) =1.499 
XU 141 =150 
C 

mil (6, •) mm mirim omsirr 

Rim, •urn 

mine, •! XNIIR MtONABU S1R1SS ' 

Rim, •ISIGAll 

mn (6. •) mm disirid difucuon' 

Rim. • I DIF All 

BASI=V0 

HGHMO 

NX--50 

mil (6, '! INIIR APPl IID LOAD ' 

Rim, • if 


mins, •) xniir noouius of iiAsncur 

Rim, • mi 

XINRUBASIUMDl/ll. 

II)[I=P , IXD ,, S)/(S , XH0U , X1NR1I 

SIRCHX- 0.0 

DO 10 1=1 NX 

sm i)=o.o 

io mm 
c 

IPRlNhS 

numu 

INF 0=0 

11=0 

100 ( 'All DOllMO.nim. IPRlNl NUV.NCON.X XI, XU, 

i obj. mm. g. rprh m, ux. m mx. nr ha r/ 

IF (XII F 10. XM) XtlhXIll-. 0000, \ 

IF MS). 10 XI4II X(SI=X(4)-. 00001 
iiuin 

IF I INFO 10.0)601020 
XI=(X(I/-X/2)I/XD 

X 2 =X( 2 ) 

XS=(X(SI-X(4)I/XD 

X 4 -M 41 

mil (8 85)11. xm XM.X(SI.X(4I 
85 f (Ml 114 41 15. 6/ 

0BJ= (RHO’XD/S) •(Xi4) ’X (2) <X( I ) ’X(SI) 
h(RH0 , X016) , lX(4)‘X(lhX(S) , XI!)l 
DO 50 J= I. NX 

SIRISS=6 , P , (XD-XII/IIXI , XUX?) , I!XS , X1*X4)* , 7)) 

smui-siRiss 

[ 

CALL BfFllC(XI,X!,X5.X4,ADP Dill. SIOPl.XI.XfW.P.XD 1 
C 

IF (XI. 10. XDi WIN 

IlFFtlAX=DFft 
IF UDIF. GlDIFliUilN 
DO 60 XU. NX 

XSMXtSh 07) X(4))/XD 

SIRlSUOflXD-XD/UXI •XhXli •( (XS*XbX4) ••!)) 

CAIl DFFIICIXI XI XS. X4. MSP Dill l, SIP), XI. XOOO P XDI 

X(SM(XIS)-.07l-X(4))/XD 

SIR! Sl-H’PUXUXl UUX) •.<!> X?) •UXS , XbX4) ••!)) 

CAll 01 Fill (XI. XI, XS. X4. MP. DIFll, Si PI, XI. XOOO, P. XD) 
BXFL=iOfFl nOFFLD/l 

if (xi io xdi Dffmm 

SlOPl-iSLPhSLPD/l 

S!RFSS=tSTR(SbSlRfS2)/2 

S!Rtl(J)=SIR!SS 
IF (XI. 10X0) MN 

GUUDIFHDIFMI '/ 001)1/(011 Ml •! 001) 

GISH (Dlf All •. 999) DIFII/IDIFAU 999) 

iNDir 

IFISIR1SS 1 1. SJRCiX I MN 
11 INCH 101)601061 

smxsmsi 

6(2) ■ ISIRflAX-SIGAl l I/S16MI 


Table 1 . Fortran program 


SPOPT (continued) 


NCH I 

mu 

SI XMI-XDiNX 

MHttm 
60 mm 

GO JOTS 

mu 

mu 

ujxi.fo.xouhfn 

mmm 

G/IMDFFHDFFALI '/. OOIIJ/IDFFAU •!. 0011 

cvMimi ' mmi/imi •. m 

mu 

umssAisnmnnFo 
if mm. I m jo 6? 
sum=sjnfss 

G/pu/sj/m-smu/siGAii 

mu 

mu 

simi-wm 

sum-mss 

socauiHUf 

c 

c ••moon a sau jo find moon if sirfss js strict t r 

C INCRFASIN6 AT JHF FIXFP FNO OF JHF LFAF SPUING "" 

C 

um.FQ.oum 

l AS JUT- 1 
00301=1.1 AST 

somsrmi 

JOIN 
JF1RSJ--IU 
00 40 J=JF INST. NX 

UlSfUN. IF. STR/UJUGO JO 40 

SHIN=STMIJ) 

JHIN=J 
40 CONTINUF 

SJun/jniNi-sjmn 

sjmiusnix 

30 CONTINUF 

sjm-STmx) 

ci?i=isjm-smu/smi 

FNDIF 

C 

C 

00 45 1=1. NX 

: until i i =o.o 

45 CONTINUF 
75 XI=X0 

srm--o.o 

NCHX=0 

mu 19. •nJMFfiAx.sjmx 

GO JO I 00 

TO STOP 

FNO 


Table 2. Fortran program SPTRIAL. 


mien rnn mam ia-h. o-d 

DIHERSIOR ADP/4.5) 

cptHtwwFim 'C: mwwmM-i 

ICOURTS 

tam=o 

HXU 5 

mm, •) erter allorable stress' 

Rim, vsigall 

mTE/6. V 'ERTER HODULUS OF ELASTICITY' 

REAIKS, VXHOD 

mTE/6. V ’ERTER APPLIED LOAD' 

READ IS, VP 

miElS.V 'ERTER BEAU LER6TH' 

READ/S, VXL 

mTE/6, V 'ERTER DERSITY' 

READ IS. VRHO 

mTE/6 V 'ERTER DESIRED HAXIHUH DEFLECT ICR' 

READ IS, VEDEF 

TEHPSML *XL •XL)/I5 , XH0D , EDEEI 

mTE/6. V 'ERTER IRITIAL GUESS FOR THICt/RESS RATIO' 

READ/S. VXRAT 

mET/mi 

TEno-.com 

mTE/6. V 'ERTER LERCTH TO VW1H RATIO' 

READ/S, VRR 

R0=XURR 

mm. am 

BASE-RO 

HGHT-TO 

miE/8, 200)BASE,HGHT 

ROD FEME //IX, 'THE IRITIAL GUESS FOR THE CROSS SECTIORAL AREA './'AT 
/THE BASF /Ti/mST PARTI a THE TAPERED BE AH IS: '. / ' RIDTH * EG 
25, ' t IHKXRESS = '.E6.5I 
XIHRTSASE‘m"5l/L2. 
lDEhWXL H 5)//5 , XHOO , XIRRT) 
mTE/8. 1001 IDEE 

100 EEMT/r THE DEELECTIDR FOR A C0HSTAR1 X-SECTION BEAD',/' RITH TH 
I ESI CIDERS WHS IS '.E6.5I 

TSTWKLVOmiW) 
miE/8, 10IITSTR 

101 EORHAT//' THE HAX. STRESS FOR A CORSTART R-SECTIOR BLAH'./' RITH T 
IHESE DIHEHS IORS IS ‘,110.51 

m. TE 18.2021 

202 EORHAT//' THE FOLLORIRG RESULTS RILL TRDICATE THE COHBIHATIOHS a 
IRE. iO. TE. ARD TO ',/' THAT PRODUCE THE DESIRED DEFLECT IOR RITHIH 
?*/- SI. AT THE ERD OF'./' EACH RUR, HESSAGE STATEHERTS HIGHT APPf 
5AR RECOttlERDIRG mflCAUCRS',/' TO A PREVIOUS RUR IR ORDER TO IR 
4CREASE ACCURACY. THE SLOPE VALUES ARE GIVER'./' IR RADIARS ARD IR 
SDICATE THE SLOPE OF THE BEAD AT THE TIP RHERE THE',/' DEFLECT IOR I 
6S GREATEST. 7! 

TEIR=TE 

D04SH- 1,4 
RRAhRE/RO 
THEIR 

DO 49 1 -- 1.80 


RMS 

st nemo 

XML 

DO SO J = i'.RX*l 
XI I /-IRE RDI/XL 
XI2PRD 

X/5MTE-T0I/XL 

X/4.I-T0 

SIRESS-6'P' IXL - Yli/imn •XUX/2H •IIX/5) W/W "21! 

c 

CALL DEE/X.AIf Dfh . SLOPE. XI. XHOU.P.XLI 
C 

IFIXI (O.Xl/TH/R 
IF/TDEE.GT DEEL/THER 
00 60 X-I.RX*/ 

XL 51- // Tt 07/ -TOl /XL STRESIS’PVXL - 
XIi/<iX/li , Xh/!2))V/X/II , XUXi4))”2ll 
CALL DEE/X. Alf HEEL /, SLPI. XI. XHOU, P. XL I 
X/5/-/ITE- 02H0I/XL STRES7S*P*/XL- 
XI)/IIXll) *XM 121 IV IX 151 •XM/4))“2)I 
CALL DEE/X ADP. HEEL 2. SLP2. XI XHOD. F, XL I 
DEEL : IDEEU*DEfL7U2 
SLOPE-- /SLPUSLP7U2 

STRESSUSIRFSISTRES2I/2 
IF/XI.EQ.XIlIHtH 
DEIH/X--DEFI 
SUM SLOPE 

RATE -P/DEEHAX 

ERDIE 

C 

IF /STRESS. L T. STRCHX I THER 
IE/RCHX E0. 1 ICO TO 61 
STRHAX- STRCHE 
RCHN 

XHt-XbXLiRX 

ERDIE 

C 

61 XI -XI -XL /RX 
SUM-STRESS 

IEIX.EO.RXVITHER 
IE/RCHX. RE. IITHER 
SUM-STRESS 

XHE-XI 

IFIXI. IT. O.OIXHIS.O 

ERDIE 

ERDIE 

60 CCH TI RUE 
GO TO 20 

ERDIE 

ERDIE 

C 

IE /XI. (0. XL) MR 
SUM-SLOPE 

DEfm-DEEL 
RATES /DEE HAX 
ERDIE 


Table 2. Fortran program SPTRIAL (continued) 


msmsmsmtim 
iFinoK.to.ua) me? 
sam-sm 

XIH-XML/XX 

HCHE-I 

mu 

s/xi-xi-xim 

sirchx--siress 

IfU.EO.XXWIHEX 
Um XE 1 1 MX 

ummss 
u I xi. u. o.o) mo 
mu 
mu 

oo mum 

n mu um xeicht- iwu / > - •a m> ■> v • ' m*xim 'now > it w 

yuf-mmm-mi) 

If I 01 fflf.il MX 

i cm i- 1 

mmumsKMi 

mms.m 

jutmmtimtmtt •) 

mis. mi 

m frnur offiicm xt xo n io 
i urns siopf v 
mints. /Mm. he, do. u. io. sum. mo 
mu is, tos weight 

uk frnur iHf xrm of m dew win mu oimioxs is •.(?. w 

intit iiniinuiinmttnnininmnnintnininitiitHiitiiin 

/? fomwx.7ifios.ixii 

in sum. si. sum urns 

XOUXJ I 
CO 10 16 

mu 

mu 

// fommur/.s, ixu 

THf-UIX/80 
49 COXIIXUf 
4 6 m-M-./S 

mmm 

48 COXIlXOf 
IKICOan.fQ.OiMX 

mu 18 . 70 /) 

m trniu mmm m w.mt no mi sue ,,nn 

!••••• WIDE WAS IRC IHf MUM THitMSS RATIO "'"'7 
CO TO SCO 

mu 

if mm co. ijmx 

mu IS. 70 !) 

70i fomn /'•••••• moo Aii stresses exceed the auouabie """ 

//'""" m IXCREASIXC THE IXIIIAl IHiaXfSS RATIO """7 

EXOIF 

TOO SIDE 

EXD 


sequence approach. The importance of having adequate initial values of the design variables will 
also be demonstrated. This will be done by comparing the stress, deflection, and weight of initial 
geometries with the optimized configurations. It will be shown how small deviations from a 
feasible initial design can result in nonconvergence of the optimization problem. 

The problem to be solved is as follows: 

• Design the minimum weight leaf springs of a suspension system that will provide protec- 
tion to the mass of 3.200 lb at a frequency of 2.2 Hz with a factor of safety of 1.4. The 
maximum G-load (load magnification) that will occur during the ascent and/or descent 
mission is 2.63. 

The first step is to determine the required spring rate for each spring (cantilever beam). From 
equation (3). one finds 


k, = f = (2.2)~4tt^ — 1.382.41 lb in (44) 

Since the isolation system has been designed such that the springs under design are all in parallel, 
one can obtain the spring rate for each beam by dividing the total spring rate by the total number 
of beams (in this case 4). Thus, the spring rate for each beam is 


A, 


A/ 

/v,„ 


1,582.41 

4 


395.6 lb/in 


(45) 


The next step is to determine the maximum load per beam. This will be done by including the 
G-load and factor of safety in the calculations. 


W ssh G(FOS) = (3.200)(2.63)( 1 .4 ) = ? ^ ^ 
N xp 4 


(46) 


where 

Wys/r = the total weight to be isolated 

G = the G-load 

FOS = the factor of safety 

N sp = the number of springs. 


Finally, the maximum deflection of the beam under the design load can be calculated from the 
definition of spring rate of a beam 
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(47) 


8 


max 


k. 


2.946 

395.6 


7.45 in 


I he beam dimensions obtained during the preliminary design phase are found in figure 4. 
Figures 5 and 6 show del lection and stress plots of this design. Inspection of the data shows that 
the maximum deflection ol this beam with the applied load is 6.75 in and the maximum bending 
stress is 102.630 psi. II the material selected for the beam is titanium Ti-6AI-4V (allowable bend- 
ing stress 104.000 psi and Young's modulus of 16E6). one can see that the stress constraint is 
met. 1 he deflection, however, yields a spring rate of 436 lb per inch, a difference of 10.2 percent 
I rom the desired 395.6 lb per inch. This difference results in a natural frequency of 2.31 Hz. a dif- 
ference ol 5.0 percent I rom the desired 2.20 Hz. This error is equal to the goal of 5.00 percent 
allowable variation I rom the design, therefore, the preliminary design configuration was deemed 
acceptable. 

The author will now proceed to design the beam using the SPTRIAL and SPOPT programs. 

Table 3 shows the input information required by SPTRIAL. 

SPTRIAL is a program to obtain initial feasible designs. It does, however, require that the 
user have knowledge of the effects of changing certain variables. For example, the thickness ratio 
has a greater effect on the maximum stress than the length-to-width ratio for a given deflection. 

This means that if a design is close to a desired deflection but the stresses are slightly above the 
allowable, it is recommended to change the length-to-width ratio (instead of the thickness ratio) to 
modify stresses without significantly affecting the stiffness of the beam. The program aims toward 
a desired deflection by varying the initial thickness ratio. The stresses for several width ratios 

are printed along with the weight of the beam and comments on whether the stresses 
exceed the allowable. The author has noted, however, that generally the designs that result from 
SPTRIAL are accepted by SPOPT to yield adequate final designs which meet both the stress and 
deflection constraints. 

Table 4 shows a listing of the results from SPTRIAL using the input data from table 3. 

Table 5 is a listing of the input data required by SPOPT. The sample data shown is from 
the first initial design of table 4. 

Table 6 shows the output listing from SPOPT for a typical optimization run. In the case 
shown, the input data from table 5 were used. The majority of the output listing is from the DOT 
optimizing code with exception to the stress and deflection values for the optimized beam. 

Even though SPTRIAL greatly helps in the selection of initial values for the optimization 
process, it cannot guarantee convergence to a feasible design every time. Table 7 shows the 
optimized results obtained for various initial designs as generated by SPTRIAL. Notice that there 
were still two cases where no feasible design was obtained. This could possibly be corrected by 
modifying some internal control parameters within DOT. but. due to the fairly consistent values of 
the optimized weights, it was felt that modifications would not improve the results greatly. 
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Beam Deflection 
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Figure 4. Initial leal spring design configuration (half length). 
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Bending Stress 



Figure 6. Stress distribution for initial leal spring design. 
Table 3. Sample of input data required by SPTR1AL. 


Enter allowable stress 

104.000 psi 

Enter Young's modulus 

I6E6 psi 

Enter applied load 

2.946 lb 

Enter beam length 

29.25 in 

Enter beam density 

0.16 Ibm in 

Enter maximum deflection 

7.45 in 

Enter initial guess for thickness ratio 

0.55* 

Enter lenath-to-width ratio 

4.5* 


* The thickness ratio is the main variable used in SPTRIAL to 
obtain the desired spring rate. Values between 0.5 and 0.75 are 
recommended as first guesses. SPTRIAL will vary this quantity 
as necessary to obtain the desired deflection. 

During the derivation of the Euler-Bernoulli beam equation, 
a major assumption is that the beam have a fairly slender 
geometry 1 1 1. The recommended minimum length-to-width ratio 
is 4.5 where the width is taken as the average between the 
widths at the free and fixed locations. 
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Table 4. Output listing from SPTRIAL. 


the initial ms for the cross see i tonal area 

HI THE me ITHimT PARI) Of THE TAPERED BEAD IS: 

VIDTH -- 6. 500 A memos -- . 905 

m omccuoH for a mm x-seciicn bead 

VITH THESE DIDENSIONS IS 5.828 

THE m. SIBFSS FOR A CONSTANT X-SECIICN BEAT! 

RUB IRESF OMRS IONS IS 92159.852 

JHF FOL LOVING RFSULJS URL INOICAIF THE CODBINATIONS Of VI. VO IF. ABB JO 
IBAI PROOOCF IHF DFSIRFO DFFIFCIION VIIBIN */- 51 AI JHF END Of 
EACH RON. HFSSAGF STATEflENTS HIOHJ APPEAR RECLWENDING MFICAJIONS 
JO A PREVIOUS RUN IN ORDER JO INCREASE ACCURACY. JHF SLOPE VALUES ARE GIVEN 
IN RADIANS AND INDICATE THE SLOPE Of THE BEAD AI THE IIP VHIRE JHL 
DEFLECTION IS GREATEST. 




DEFLECTION VE VO TE TO STRESS SLOPE 
. 259486*0/ . GSOOOEm . 650001*01 . 550616*00 . 904806*00 . 1022/6*06 . 468921*00 

THE HEIGHT Of THE BEAD VI TH THESE DIDENS IONS IS 19. 095 




DEFLECTION VE VO TE TO STRESS SLOPE 
. 75/5 1 6*01 . 48/506*01 . 650006*01 . 40/ L 66*00 . 904801*00 . 10647E*06 . 461/51*00 

THE HEIGHT OE THE BEAD VI TH THESE DIDE NS IONS IS 1/ 2/6 

uminimiunmnumtnuntminmu»tuimtmumunnnimni 




DEFLECTION VE VO TE TO STRESS SLOPE 
. 258006*01 . 52500E*0I . 65000E*0L . 475026*00 . 904806*00 . 111576*06 . 4600/6*00 

THE HEIGHT Of THE BEAD VI TH THESE DWENSIONS IS 16. 285 




DEFLECTION VE VO TE TO STRESS SLOPE 
. 75674E*01 . 162506*01 . 650006*01 . 565506*00 . 904801*00 . 118086*06 . 464256*00 

THE HEIGHT Of THE BEAD VITH THESE DIDENSIONS IS 14. 622 


****** SODE OR ALL STRESSES EXCEED THE ALLOVABLE ****** 
****** TRY INCREASING THE INITIAL THICXNESS RATIO ****** 


Table 5. Sample of input data required by SPOPT. 


Enter the beam length 

29.25 in 

Enter initial guess for uy 

6.50 in* 

Enter initial guess for uy 

6.50 in 

Enter initial guess for 

0.3506 in 

Enter initial guess for /„ 

0.9048 in 

Enter material density 

0. 1 6 Ibm in ' 

Enter allowable stress 

104.000. psi 

Enter desired deflection 

7.45 in 

Enter applied load 

2.946 lb 

Enter Young's modulus 

I6E6 psi 


In order to preclude computational difficulties arising from 
design cases where uy is equal to uy. and when t, is equal to 
SPOPT adjusts the input information to eliminate the possi- 
bility of a singularity. For this case, uy is set equal to 6.4999 
without significantly affecting accuracy. 


It should be noticed from table 7 that all the initial designs meet the deflection constraint of 
7.45 in within approximately 0.050. The stress constraints, however, are violated many times but 
this does not preclude convergence to a relative optimum. This indicates that the initial designs 
need not be feasible in order for the problem to converge. All final designs were within 0.015 in 
of the desired deflection, and the stresses were within 0.5 percent of the desired stress. 

Figures 7 through 10 show the convergence history of the optimum design variables, 
optimum deflection, optimum stress, and minimum weight, respectively, for the initial design beam 
of figure 4. Notice that the final weight is approximately I lb heavier than the initial design (table 
7 ). Increases to the base dimensions uy and /„ were made by DOT in order to obtain a solution 
closer to the constraints without significantly violating them. 

It is interesting to find out what the optimum configuration would be if manufacturing (i.e.. 
material availability) or allocated space restrictions were included in the optimization routine. In 
figure I I it has been assumed that the only titanium available with the desired properties is a plate 
with a thickness of 0.930 in. If 0.020 in is allowed for machining, this means that the maximum 
material thickness available is 0.910. It has also been assumed that, due to space restrictions (i.e.. 
to prevent interference with adjacent hardware), the maximum width allowable is 6.50 in. These 
limitations are very close to actual restrictions during the preliminary design effort and limit the 
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Table 6. Output listing from SPOPT (DOT optimizer). 


om om rmm integer paraoeters 


D 

D 

0 0 

T 

mm -- 

0 

61 NCOLA -- 

8 

in IPRNTl -- 

0 

D 

D 

= o $ o ---- 

T 

21 ISCAL -- 

4 

7! KM -- 

0 

121 IPRNT2 -- 

0 

D 

D 

0 0 

T 

5i im -- 

40 

81 JTOAX -- 

20 

151 MITE -- 

0 

om 

om 

T 

41 1 TROOP - 

2 

91 1 TRUST = 

2 







51 MITE -- 

7 

101 JPRINT - 

0 




DESIGN OPTIOI2ATION MS 


1C) COPMGHI, 1985-89 
ENGINEERING DESIGN OPE MW 108, INC 


storage reouireoents 

ARRAY DWENSION REQUIRED 
H 800 20? 

m 200 si 


AIL RIGHTS RESERVED, NORLDVIDE 


- INITIAL VARIABLES AND BOUNDS 


VERSION 2. 00 LOVER BOUNDS ON THE DECISION VARIABLES (XL-VECTOR! 

II 9.99000E-01 1.000001*00 2.49900E-0I 2.50000E-0I 

DECISION VARIABLES (X-VECTORI 

- YOUR INTEGRITY IS OUR COPYPROTECTION - I) 0.499991*00 0500001*00 5.5001 OE-OI 9.04800E-0I 

UPPER BOUNDS ON THE DECISION VARIABLES IXU-VECTORL 
II 9.99900E*00 l.00000E*0l l. 49900 E *00 l.5mE*00 


CONTROL PARAOETERS 


- INITIAL FUNCTION VALUES 


OPTTOI7ATION OEM. 

OEM -- 

1 

NUOBER Of DECISION VARIABLES, 

NDV-- 

4 

RUBBER OF CONSTRAINTS, 

NC0N-- 

5 

PRINT CONTROL PARAOETER, 

IPRINT-- 

2 

GRADIENT PARAOETER, 

h TRAD-- 

0 


GRADIENTS ARE CALCULATED BY DOT 
XLIH .999 XUTII 9.999 
XU?) 1.000 XU 12! 10.000 
XL! 51 .250 XU 1 51 1.499 
XU 41 .250 XUT4I 1.500 


OBJ = 19.095 

CONSTRAINT VALUES 16-VECTOR! 

1 1 -8.59I52E-05 - 1.6585 IE-02 040629E-05 

- BEGIN CONSTRAINED OPT U1I 2 AT ION: NED TIE THUD 

-ITERATION I OBJ - /. S8950E*0I 


- SCALAR PROTM PARAMETERS 


DECISION VARIABLES (X-VECTORI 
1 1 6.42979E*00 6.56565E*00 5.64799E-0I 8.98225E-0I 


REAL PARAOETERS 
II CT =-5.tmE-02 
2ICTTIIN -- 5.LXJ000E-05 
5) DABOBJ - I.90948E-02 
41 DELOBJ - I.00000E-05 
51 DOB J! = I.00000E-0I 
6ID0BJ2 -* 5.818951*00 
71 DXI -- I.00000E-02 


8J DX2 -- l.50000E*00 
91 FDCH -- I.00000E-05 
lOIEDCffl -- I.00000E-04 
III RHVLH2 -- 4.000001-0! 
121 DABS TR: O.OmE*DO 
151 DELSTR = I.00000E-05 


-ITERATION 2 OBJ - l.85597E*0l 

DECISION VARIABLES (X-VECTORI 
1 1 6. 2 I 285 E *00 C.I7I97E*00 5.59267E-0I 9.28420E-0 1 
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fable 6. Output listing from SPOPT (DOT optimizer) (continued) 


-imriOH 5 OBJ -- l.85597E*0l 
decisior mimes ix-vectori 

1 1 6.212856*00 6.(7(976*00 5.59267E-0I 9.284206-0! 


FURCTIOR CALLS - 46 

******* THE DEELECTIOR FOB THE OPT MM! BEAT! IS ******* 7. 46416 
******* THE STBESS FOR THE OPTIMO BEAU IS ******* 104429. 


- ITERATION 4 OBJ -- 1.855976*01 

DECISIOR MIMES (HECTOR) 

I) 6. mm 6. 171971*00 5. 59267 E -01 9.28420E-0I 


- OPTHimTICH IS COMPLETE 

RUBBER OF ITERATIORS - 4 

COHSTRAIHT TOLERARCE. CT = -5.00000E-05 

THERE ARE 5 ACTIVE CORSTRAMS ARO 0 MATEO CORSTRAIRTS 
COHSTRAIHT ROBBERS 

1 2 5 

THERE ARE 0 ACTIVE SIDE CORSTRAIRTS 
TERHIHATIOH CRITERIA 

RELATIVE C0RVER6ERCE CRITERIOR HAS LIFT FOR 2 COHSECOTIVE ITERATIORS 
ABSOLUTE CORVERBERCE CRITERIOR HAS BET FOR 2 COHSECOTIVE ITERATIORS 

- 0PTM7ATI0R RESULTS 
OBJECTIVE, FIX) -- l.ffifl 

DECISIOR VARIMES, X 

ID XL X XU 

/ 9990001-0! 6. 2(285F*00 9.999001*00 

2 I.00000EHX) 6.171971*00 l.00000E*0l 

5 2.49900E-0I 559267E-0I (.499001*00 

4 2.600006-0! 9.28420E-0I (.600006*00 

CORSTRAIRTS, 0(XL 

I) 7. 91 28 IE-04 U6577E-05 -2.79490E-05 
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Design Variables 


Table 7. Optimization results for various initial designs. 


Design 

\v\. 



t„ 

Weight 

Del 

Stress 

Iteration 

Initial 

6.500 

6.500 

0.351 

0.905 

19.095 

7.395 

102.270 


Final 

6.213 

6.172 

0.339 

0.928 

18.360 

7.464 

104.429 

46 

Initial 

4.875 

6.500 

0.407 

0.905 

17.776 

7.375 

106.470 


Final 

5.193 

7.002 

0.402 

0.878 

18.593 

7.435 

104.516 

66 

Initial 

3.250 

6.5(H) 

0.475 

0.905 

16.285 

7.380 

1 1 1 .570 


Final 

3.520 

8.122 

0.493 

0.817 

18.422 

7.459 

104.022 

71 

Initial 

4.000 

6.500 

0.490 

0.910 

17.609 

6.750 

102,630 


Final 

4.245 

7.730 

0.467 

0.829 

18.662 

7.460 

104.470 

109 

Initial 

5.850 

5.850 

0.286 

0.995 

17.541 

7.427 

108.890 


Final 


***** No Feasible Design Was Obtained 

***** 


Initial 

4.388 

5.850 

0.336 

0.995 

16.321 

7.472 

1 15.660 


Final 

5.294 

6.977 

0.390 

0.884 

18.610 

7.436 

104.472 

88 

Initial 

2.925 

5.850 

0.41 1 

0.995 

15.100 

7.368 

121.600 


Final 


***** No Feasible Design Was Obtained 

***** 


Initial 

1 .463 

5.850 

0.498 

0.995 

13.624 

7.407 

133.830 


Final 

1.516 

8.799 

0.625 

0.767 

17.203 

7.438 

104.517 

111 



Figure 7. Flistory of design variables versus number of iterations. 
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Deflection Constraint 



Number of Iterations 

Figure 8. History of deflection constraint versus number of iterations. 



Number of Iterations 

Figure 9. History of stress constraint versus number of iterations. 
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Number of Iterations 


Figure 10. History ot weight versus number of iterations. 
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Figuie I I . Design space tor leal spring with additional side constraints. 
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number of potential optimum solutions. However, it can be seen that the weight is essential l\ the 
same as the optimum solution from table 7 (18.662 lb), and a good improvement toward meeting 
the desired deflection and stress constraints is obtained. The minimum weight under these con- 
ditions is 18.664 lb with the follow ing parameters: 

ny - 6.09 in 
ic,, - 6.50 in 
-0.351 in 
t„ - 0.910 in 
Deflection - 7.46 in 
Stress - 103.950 psi. 


TEST RESULTS 


Once the preliminary design was completed and a geometry selected, a spring was manufac- 
tured from 6061-T6 aluminum alloy to verify the configuration. Since the final configuration was 
to be manufactured out of an expensive titanium alloy, the decision to proceed would be based on 
the outcome of this test. 

The test parameters were as follows: 

• Maximum load per cantilever - 800 lb 

• Maximum expected deflection - 2.934 in 

• Maximum expected stress - 27.868 psi 

• Expected spring rate - 273 Ib/in. 

Strain gauge and displacement indicator locations for the test hardware were as indicated in figure 
12. Test procedures and results are in references 8 and 9. Pertinent information is summarized 
below: 


• Maximum applied load - 800 lb 

• Maximum measured deflection - 3.052 in 

• Maximum measured stress - 27.160 psi* 

• Measured spring rate - 262 lb/in. 

It should be noted that the maximum measured stress was obtained at gauge number 4 in figure 
12 (9.5 in from the clamped edge). The actual calculated location of maximum stress is at 10.38 in 
from the clamped edge. The calculated stress at 9.5 in is 27.850 psi. a difference of 2.47 percent. 
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INSTRUMENTATION LOCATIONS 
TENSION SIDE 

O EDI LOCATIONS 

□ UNIAXIAL STRAIN GAGE LOCATIONS 



HALF OF Z-SPRING 
TENSION SIDE 


FIGURE 12 


Figure 12. Location of deflection and strain gauges. 
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Figure 13 shows the predicted versus the test values of the deflection for the aluminum test 

beam. 


CONCLUSIONS 


The results of this study show that although time consuming trial-and-error iterations were 
performed during the initial design of the leaf springs for the HST/SSE. the resulting design was 
very near an optimum design for the configuration analyzed. The study also shows that with the 
availability of personal computer-based optimization software, fairly complicated problems can be 
handled with fast solutions and reliable final designs. It is interesting to point out that constraints 
and limitations such as material availability and possible interference with adjacent hardware can be 
included in the optimization procedure as mathematical constraints on the numerical minimization 
problem. 

It is important to note that although the derivation of the deflection equations for the spring 
was based on small deflections, comparison with nonlinear finite element solutions show that for 
the range of deflections required, the difference between both solutions is acceptable. Care must be 
exercised in order to justify the linear approximation for applications with larger deflections. 

The time has come for design engineers to take advantage of the powerful tools available 
for developing lightweight and structurally sound hardware. All that is required is the desire to 
learn and the awareness that the state-of-the-art is advanced by inquisitive minds. 
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Beam Deflection 



Applied Load 


Figure 13. Test results for leaf spring (6061-T6 aluminum). 
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